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We study the third order solutions of the cosmological density perturbations in the Horndeski’s 
most general scalar-tensor theory under the condition that the Vainshtein mechanism is at work. 
In this work, we thoroughly investigate the independence property of the functions describing the 
nonlinear mode-couplings, which is also useful for models within the general relativity. Then, we 
find that the solutions of the density contrast and the velocity divergence up to the third order ones 
are characterized by 6 parameters. Furthermore, the 1-loop order power spectra obtained with the 
third order solutions are described by 4 parameters. We exemplify the behavior of the 1-loop order 
power spectra assuming the kinetic gravity braiding model, which demonstrates that the effect of 
the modified gravity appears more significantly in the power spectrum of the velocity divergence 
than the density contrast. 

PACS numbers: 04.50.Kd,95.36.+x,98.80.-k 


I. INTRODUCTION 

The accelerated expansion of the Universe is one of the most fundamental problems in modern cosmol¬ 
ogy. The standard cosmological model introducing the cosmological constant is consistent with various 
observations [1, 2], However, the small value of the cosmological constant raises the problem of fine tuning 
[3-5]. As an alternative to the cosmological constant, the cosmic accelerated expansion might be explained 
by modifying gravity theory, e.g., [6-18]. In the present paper, we focus on the most general scalar-tensor 
theory with the second order differential field equations [19, 20], which was first discovered by Horndeski 
[21]. The Horndeski’s most general scalar-tensor theory, including 4 arbitrary functions of the scalar field 
and kinetic term, reduces to various modified gravity models by choosing the specific 4 functions. Because 
the Horndeski’s theory includes a wide class of modified gravity models, we adopt it as an effective theory 
of the generalized theories of gravity. 

In the present paper, we investigate the aspects of the quasi-nonlinear evolution of the cosmological 
density perturbations in the Horndeski’s most general scalar-tensor theory assuming that the Vainshtein 
mechanism is at work [22-25]. The Vainshtein mechanism is the screening mechanisms, which is useful to 
evade the constraints from the gravity tests in the solar system. We investigate the effects of the nonlinear 
terms in the matter’s fluid equations as well as the nonlinear derivative interaction terms in the scalar 
field equation. In a previous work [26] , the second order solution of the cosmological density perturbations 
is obtained. In the present paper, we extend the analysis to the third order solution, which enables us to 
compute the 1-loop order matter power spectrum. 

There are many works on the higher order cosmological density perturbations and the quasi-nonlinear 
matter power spectrum, which have been developed from the standard perturbative approach (see e.g., 
[27-37]). Improvements to include the non-perturbative effects have been investigated, e.g., [38-43], but 
we here adopt the standard perturbative approach of the cosmological density perturbations as a starting 
place for the analysis of the Horndeski’s most general scalar-tensor theory. Related to the work of the 
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present paper, we refer the recent work by Lee, Park and Biern [44] , in which a similar solution is obtained 
for the dark energy model within the general relativity. 

This paper is organized as follows. In section 2, we review the basic equations and the second order 
solution [26]. In section 3, we construct the third order solutions of the cosmological density perturbations. 
Here, we carefully investigate independent functions of mode-couplings describing nonlinear interactions. 
In section 4, we derive the expression of the 1-loop order power spectra of the matter density contrast 
and the velocity divergence. In section 5, an expression for the trispectrum for the density contrast is 
presented. In section 6 , we demonstrate the behavior of the 1-loop order power spectra in the kinetic 
gravity braiding model. Section 7 is devoted to summary and conclusions. In appendix A, definitions 
of the coefficients to characterize the Horndeski’s theory are summarized. In appendix B, definitions of 
the functions to describe the nonlinear mode-coupling for the third order solutions are summarized. In 
appendix C, a derivation of the 1-loop power spectra is summarized. Expressions in appendix D are useful 
for the deviation of the 1-loop power spectra. Appendix E lists the coefficients to characterize the kinetic 
gravity braiding model. 


II. REVIEW OF THE SECOND ORDER SOLUTION 


Let us start with reviewing the basic formulas [23, 26]. We consider the Horndeski’s most general 
scalar-tensor theory, whose action is given by 


S = 


d A x y /~ z g(C GG + C m ), 


( 1 ) 


where we define 

£gg = K{(j>, X) - G z {<f>, X)\2<j> + G 4 (& X)R + G iX [(P </>) 2 - (V M V^) 2 ] 

+G 5 (0, X)G^WV^ - i G 5X [(D0 ) 3 - 3D0(V M V,,0 ) 2 + 2(V Al V^) 3 ], ( 2 ) 

where K, G 3 , G 4 , and G 5 are arbitrary function of the scalar field c\> and the kinetic term X := — (dcj)) 2 / 2, 
Gi X denotes dGi/dX , R is the Ricci scalar, Gis the Einstein tensor, and £ m is the Lagrangian of the 
matter field, which is minimally coupled to the gravity. 

The basic equations for the cosmological density perturbations are derived in Ref. [23]. Here, we briefly 
review the method and the results (see [23] for details). This theory is discovered in [19] as a generalization 
of the galileon theory ([45], see also [20, 46-64]), but the equivalence with the Horndeski’s theory [21] 
is shown in [20]. We consider a spatially flat expanding universe and the metric perturbations in the 
Newtonian gauge, whose line element is written as 

ds 2 = —(1 + 2 $(t,x))dt 2 + a 2 (t)( 1 — 24/(f, x))dx 2 . (3) 

We define the scalar field with perturbations by 

0 +5(j>(t,x), (4) 


and we introduce Q = H5<j)/(j). 

The basic equations of the gravitational and scalar fields are derived on the basis of the quasi-static 
approximation of the subhorizon scales [23] . In the models that the Vainshtein mechanism may work, the 
basic equations can be found by keeping the leading terms schematically written as (ddY ) n , with n > 1, 
where d denotes a spatial derivative and Y does any of <5, 4' or Q. Such terms make a leading contribution 
of the order (L^ddY) n , where Ljj is a typical horizon length scale, and we have 

= 2§TpO (,] + Jfp (V 2 4>V 3 Q - 3,3,*8*8*0), 


V 2 {Rt' i - Gt® - AiQ) 


( 5 ) 




and 
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StV 2 * = - *V’(J - 

- a,a,iaVQ) - gA_e |s >, 

where p m is the background matter density and <5 is the matter density contrast, and we define 

Q (2) := (V 2 <3 ) 2 - (didjQ) 2 , 

Q (3) := (V 2 Q) 3 - 3V 2 Q (Si^Q) 2 + 2 (d^Q) 3 . 

The equation of the scalar field perturbation is 


A 0 V 2 Q - A jV 2 '!' - A 2 V 2 <I> 

b 2 


Bn 


a 2 H 2 


q( 2 ) - 


a 2 H 2 

Co 

(Ah 4 


(v 2 $v 2 g - did j &&& i Q) - 


a 2 H 2 


eP, -^ <3> = »- 


(v 2 <Fv 2 g - c^cwg) 
(v 2 $v 2 t - d^TcFcT'F) 


where we define 


( 6 ) 

( 7 ) 

( 8 ) 


( 9 ) 


U (3) g(2) V 2 $ _ 2V 2 Qd l d j Qd l d j $ + 28 l 8 j Qd :i 8 k Qd k d i ^>. (10) 

Here the coefficients Tt, Ai, B i, C\, etc., are defined in Appendix A. A i: B i: and Ci are the coefficients 
of the linear, quadratic and cubic terms of T, <f>, and Q , respectively. 

From the continuity equation and the Euler equation for the matter fluid, we have the following equations 
for the density contrast 6 and the velocity field u l , 

^ 75 ^ + ~3i[(l + 5(t,x))u%x)] = 0 , ( 11 ) 

at a 

(t Cl 1 1 

— ’ + -w*(t,x) + vP (t , x) dj u 1 (t , x) = --<%$(£,x), (12) 

ot a a a 

respectively. The properties of the gravity sector is influenced through <f> in (12), where $ is determined 
by Eqs. (5), ( 6 ) and (9). 

Now introducing the scalar function 9 = Vu /(aH), which we call velocity divergence, and we perform 


the Fourier expansions for S and 0 , 

S(t,x) = (2^)3 / d3 P 5 ^ P) e * P ' X ’ ( 13 ) 

u J (*> x ) = oTl 3 [ d 3 p-^-aH9{t, p)e lp x . (14) 

(2^ J p 2 

In the similar way, we perform the Fourier expansions for $, 'F, and Q. Then, the equations for the 
gravity (5) and ( 6 ) lead to 

~p 2 (JV'i'(f,p) - g T $(t, p) - Aig(f,p)) = ^^r[t,p ; g,g] + .J^T[t,p;Q,®], (15) 

„ a 2 

-p (G T ^(t, p) + A 2 Q(t, p)) - —p m 6(t, p) 

= - 2 !ah 2 P: Q’ Q]~ - 3^4 3 i[* ,p;g»g,g]> ( 16 ) 
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respectively, where we define 


T[t,p-,Z 1 ,Z 2 ] = J dk 1 dk 2( 5®(k 1 +k 2 -p)(fc 1 2 fc 2 2 -(k 1 -k 2 ) 2 )Z 1 (t,k 1 )Z 2 (f,k 2 ) ) (17) 


l 1 [t,p-,Zi,Z2,Z 3 \ — 


dkidk 2 dk 3 <5 (3 )(ki + k 2 + k 3 - p) 


(2t r) 6 

-klk 2 2 kl + 3fc 2 (k 2 • k 3 ) 2 - 2(kr • k 2 )(k 2 • k 3 )(k 3 • k 3 ) 


Zi(t, ki)Z' 2 (t, k 2 )Z 3 (f, k 3 ), (18) 


where Z x , Z 2 and Z 3 denote any of Q, $ or 4'. The equation for scalar field perturbation (9) leads to 
—p 2 (AoQ(t, p) - Ai^(t, p) - A 2 $(t, p)) 

= p; Q- Q] + r t*. p; Q> *1 + ^§2 r [*. p; Q ,+ ^§2 r [*> p; ^ 

+ ^^4Si[t,p;Q,Q,Q] + -^ I S 2 [t,p;Q,Q,$], (19) 

where we define 


^ 2 [t, p; Z Xl Z 2 , Z 3 \ 


1 

( 2 ^ 0 ® 


dk 3 dk 2 (ik 3 (5^ 3 '(ki + k 2 + k 3 


P) 


-fc 2 /cffc 2 + (k 3 • k 2 ) 2 kl 


+ 2fc 2 (k 2 • k 3 ) 2 


2 (ki • k 2 )(k 2 • k 3 )(k 3 • k t ) 


Z\(t 1 ki )Z 2 (t, k 2 )Z 3 (t, k 3 ). 


( 20 ) 


The fluid equations (11) and (12) lead to 


1 dS(t. p) 
H dt 


+ 0(t,p) = 


1 

w? 


1 dfl(£,p) 
H dt 




dkirik 2 5 (3) (ki +k 2 


r$(*,P) 


a 2 H 2 
dkidk 2 ^ 3) (k! +k 2 


p)a(ki, k 2 )0(t, ki)5(t, k 2 ), 


p)/3(ki, k 2 )0(t, ki)6>(t, k 2 ), 


( 21 ) 

( 22 ) 


where we define 


a(ki,k 2 ) 


^(ki,k 2 ) 


1 + Kl • K2 

ki 

(ki ■ k 2 )|ki + k 2 | 2 
2 fc 2 fc| 


(23) 

(24) 


Note that a(k 3 , k 2 ) doesn’t have the symmetry with respect to exchange between ki and k 2 . We find 
the solution in terms of a perturbative expansion, which can be written in the form 


Y(t, P) = ^Tn(i,p), 

71—1 


(25) 


where Y denotes <5, 0, $ or Q , and Y n denotes the nth order solution of the perturbative expansion. 
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Neglecting the decaying mode solution, the linear order solution is written as [23, 65] 

M*,P) = D + (t)S L (p), 

Oi(t,p) = -D+(t)f{t)S L (p), 
a 2 H 2 

$i(*,p) =-p), 

pZ 

a 2 H 2 

^i(t,p) = -— D+(t)K*{t)6L(p), 

pz 

a 2 H 2 

Qi(^p) =-p), 


where D+ ( t ) is growth factor obeying 

d 2 D+(t) 

dt 2 


+ 2 H 


dD + {t) 

dt 


+ L(t)D + (t) = 0, 


with 


m 


(A 0 Tt — A\)p m 

2 (AgQ^ + 2A1A2QT + A^Tt) 


(26) 

(27) 

(28) 

(29) 

(30) 


(31) 


(32) 


and <5 l(p) describes the linear density perturbations, which are assumed to obey the Gaussian random 
distribution. Here we adopt the normalization for the growth factor D + (a) = a at a -C 1, and introduced 
the linear growth rate defined by f(t) = d\n D + (t)/ In a. 

The second order solution is written as (see [26] for details), 


(5 2 (t,p) = D%(t) (w a (p) - yA(f)W 7 (p)^ , (33) 

®2(t,p) = -D 2 + {t)f (W Q ( P ) - yA e (t)W 7 (p)^ , (34) 

2tt 2 

<f> 2 (t,p) =- 5 — D 2 + (i){n<s,{t)yV a {p) + A$(t)W 7 (p)), (35) 

pZ 

a 2 H 2 

^2(t, p) =- D+(t)(K*(t)W a (p) + A^(f)W 7 (p)), (36) 

pZ -r 

(fo{t,p) =-^-D 2 + (t)(K Q (t)W a (p) + \ Q (t)W 7 (p)), (37) 


where the coefficients kq , A, \g , A$, A$, and Aq, are determined by the functions in the Lagrangian 

and the Hubble parameter, whose definitions are summarized in appendix A. Here W Q (p) and W 7 (p) are 
defined as 


Wa(p) = 

(2 tt) 3 j 

1 dkidM^ki+ka- 

- p)a (s) (ki,k 2 ),5L(ki) ( 5L(k 2 ), 

(38) 

W 7 (p) = 

(2 tt) 3 j 

[ dk 1 dk 2< 5( 3 )(k x +k 2 - 

- Pb( k i, k 2 )5 L (ki)5 L (k 2 ) 

(39) 


with 


(ki, k 2 ) = 1 + 


ki • k 2 (fc 2 + k |) 
2 kfk 2 


(ki • k 2 ) 2 
k 2 k 2 ' 


(40) 


7 (ki,k 2 ) = 1 


(41) 
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where a^ s ^(ki,k 2 ) is obtained by symmetrizing a(ki,k 2 ) with respect to ki and k 2 , and 7 (ki,k 2 ) is the 
function to describe the mode-couplings for the nonlinear interaction in the gravitational field equations 
and the scalar field equation, a;W(ki, k 2 ), /3(ki,k2) and 7 (ki,k 2 ) have the symmetry with respect to 
exchange between ki and k 2 . One can easily check that the functions to describe the nonlinear mode- 
couplings, a( s )(ki,k 2 ),/3(ki,k2), and 7 (ki,k 2 ) satisfy 

/?(kr,k 2 ) = a (s) (k!,k2) - 7 (ki,k 2 ). (42) 


III. THE THIRD ORDER EQUATIONS 


In this section we consider the third order solutions. The third order solution of the cosmological 
density perturbations has been investigated in various models [27-33, 35-37, 44]. We present the third 
order solution for the Horndeski’s theory in the cosmological background. Our results are general and 
applicable to various modified gravity models. Plus our results are useful for the case of the general 
relativity because we clarify the independence property of the mode-coupling functions and the relevant 
parameters to characterize the third order solution. We start with solving the third order equations for 
gravity and scalar field 


- P 2 P) - p) - AiQ 3 (t, p)) = ^^2 r[t, p; Qi,Q 2 \ + (r[f, p; Qi, 




+r[f,p;<?2, , I , i]), 


(43) 


- V 2 (t, p) + A 2 Q 3 (t, p)) - C -^-p m S 3 {t, p) = --^yl’[t, p; Qi, Q 2 \ - B3 


H 2 


(r[t,p;Qi, 


2 


+r[*,p;Q 2 ,^i]) 


C! „ 


3a 4 H 4 


i[t, p; Qi, Qi, Qi], (44) 


-P {A 0 Q 3 (t,p) - Ar^ 3 (t,p) - A 2 4> 3 (t,p)) = --^2r[f,p;Qi,Q 2 ] + 


(r[<, p; Qi , 1 

+r[t,p;Q 2 , 4 'i]) + (r[t,p;Qi,$ 2 ] + P[t, p; Q 2 , $1]) + (r[t, p; 'Pi, $ 2 ] 

C, 


Bi 

a 2 H 2 
B 3 


+r[t, p; fk 2 ,+ q4 ^ 4 5i[f, p; Qi, Qi, Qi] + q4 j^- 4 [t, p; Qi, Qi, $i]- 

Inserting the first and the second order solutions into the above equations, we finally have 

T T ^ 3 (t, p) - Gr^ 3 {t, p) - AiQ 3 (t, p) =- —Dl(t)[(B 1 KQ + 2 B 3 k$kq) W 7Q (p) 

+ ( BikqXq + B 3 (k$A q + kqA$)) W 77 (p)^, 
q 2 a 2 H 2 ( 

g T ^ 3 {t, p) + A 2 Q 3 (t, p) + -^PmS 3 (t, p) = — — D\(t) i (B 2 k 2 q + 2 B 3 k^kq) W 7a (p) 

qW 5 ( P ) N 


(45) 


(46) 


(47) 


c 

+ ( B 2 kq\q + B 3 (k^Xq + kqA^ 2 ))W 77 (p) H — k 3 
AoQ 3 {t,p) - Ai^ 3 (t,p) - A 2 $ 3 (t,p) =- -^—D\(t)\A-2B 0 KQ + 2B 1 k^,k q 

+ 2B 2 k$kq + 2B 3 k$K\£^J W 7 a(p) + ^( — 2BqkqAq + Bi(k$Aq + kqX vp) 

B 2 (k$Aq + kqAq) -f- 7? 3 (k<J)A\[/ + A<j»))W 77 (p) + (Uq^q + Cik^kq)W^(p) j, (48) 
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where we define W 7a (p), W 77 (p) and W|(p) by Eqs. (Bl), (B2) and (B3), respectively, in appendix B. 
Then, the gravitational and the curvature potentials, and the scalar field perturbations are written as 

a 2 H 2 

$30, p) = —— (K$(t)<f 3 (t, p)+ (a^(t)yv 7a (p) + M$(t)yv 7 -y(p)+ ^(t)yv^(p))), ( 49 ) 

a 2 H 2 

$30,p) = —(k+O^O,p) + 0+0) (ct^0)w 7Q ,(p) + /z*0)w 77 (p) + i/4t0)w^(p))) , (50) 

Q 3 O,p) = («o0)^0i p) + 0+0) OgO)W 7 <«(p) + mqO)W 77 (p) + ^q0)w«(p))) > ( 51 ) 

where the coefficients <r$0) M*0), z'+O), etc., are defined in appendix A. The third order equations for 
63 0, P) an d 6*3 0,p) are, 


1 dS 3 {t. p) 
H dt 


+ 03 0, P) 


(2t r) 3 


dkidk 2 i5 (3) (ki + k 2 - p)a(ki, k 2 )(0i(t, ki)^ 2 (t, k 2 ) + 9 2 (t,k 1 )5 1 (t,k 2 )), (52) 


1 903 0, p) p 2 

H dt + Jp) 03 ^ ,P ^ _ ^2ff2 $3 ^’ P ^ 

" dk,dk 2 5^ (k, + k 2 - p)/3(ki, k 2 )0!0, ki)0 2 O, k 2 ). 


(2tt) 3 


Using the first and the second order solutions, these equations are rewritten as 
1 9<5 3 (t, p) 


H dt 


030, P) 


1 90 3 (t, p) , ( n , H \ n _ s P‘ 


H dt 


+ 2 + 772 0 3O, P) - $30, P) 


H 2 


a 2 H 2 


(53) 


= 0 + 0 )/ ^W Qafl (p) - yAW Q7 fl(p) + WaaL (p) - yA e W a7 i(p)^ , (54) 


= 2iA^0)/ 2 -W QQ (p) + -A e w a7 (p) + W 7 a(p) - -AeW 77 (p) , 


(55) 


where we introduce the functions defined by Eqs. (B7) to (B12), for which we find that the following 
relations hold, 


2W 7 «(P) = W« 7fi (p) + 2W 77 (p) - W,(P), (56) 

2W qq ( P ) = W aaR (p) + W aaL (p ), (57) 

2Wq 7 ( P ) = Wa 7 ii(p) + W« 7 l(p). (58) 

Then, Eqs. (54) and (55) reduce to 

i ? ^dt 1)1 + 03 ^’ = D +^ (^ 2W ““(P) ~ ^AW a7fl (p) - yA e W Q7i (p)^ , (59) 

^903|^p) + (^+^j 0 3 O,p) - ^$ 3 0, P) = 2D 3 + 0)/ 2 (-W q «(p) 

+ + i^)W ai R{p) + yAeW a7 L(p) + (l — -Ag^ W 77 (p) — -W^(p)^ . (60) 
















Combining these two equations, we have the third order equation for <5,3 (f, p) as 


d 2 S 3 (t,p) , 0 ^dSsitip) ( TU ^ U ^ 0 

- di? -^ - — -f- L{t)ds{t, p) — Ss3\t, p), (61) 

where we define 

S S3 (t) = D\(t)(N aa (t)W aot (p) + N aryR (t)W a ^ R (p) + N a ^ L {t)yV aiL {p) 


+ iV 77 (t)W 77 ( p) + Ns(t)Ws( P )), 

(62) 

and 


N aa {t)=6f 2 H 2 -2L, 

(63) 

N aiR (t) = -f 2 H 2 - ^f 2 H 2 A + 2 -L\ - pHX + X -H 2 a *, 

(64) 

N a ^ L (t) = —f 2 H 2 - *-f 2 H 2 A + 2 -L\ - ifHX + JV 7 , 

(65) 

A7 77 (f) = —2 f~H 2 + —f 2 H 2 X + — fHX + H 2 (<r$ + //$), 

(66) 

ATg(t) = f 2 H 2 + H 2 —<7$ + ' 

(67) 

where used Eqs. (A24) and (A28), and 


m = jf(~ 2 f H2 ~ L ~ f H2 - fH), 

(68) 

which follow from the definition of the growth rate f(t) = dlnD + /d\na and Eq. (31). 
A r a 7 i(t) is equivalent to N aiR (t), using (65) and (64), and N 7 (t) = \H 2 cr^, which is 
Eqs. (A22) and (A29). Then, we write 

We can prove that 
demonstrated from 

N aj (t ) = N a ~, R (t) = N ajL (t). 

(69) 


The general solution of Eq. (61) with (62) is 

<S 3 (t, P ) = c + (p)D + (t) + c_(p) J D_(t) +J* g ~ ( * )J+( ^~ < ^ +(f)J ~ ( * ,) ^ 3 (*^p)^, (70) 

where D+(t) and D-(t) are the growing mode solution and the decaying mode solution, satisfying equation 
(31), c+(p) and c_(p) are integral constants, and W(t) is the Wronskian defined by W(t) = D+(t)D-(t) — 
D+(t)D_{t). Since we assume that the initial density perturbations obey the Gauss distribution, we set 
c ±(p) = 0, as is done in deriving the second order solution. Then, the solution of the third order density 
perturbations is given by 


fa(t,p) 


D+{t) ^KS3{t)W aa (p) - yA 53 (t)W Q7 ij(p) - yA5 3 (t)Wa 7 L(p) 
- Jj-Mi)W 77 (p) + ^(t)Wj(p)^, 


(71) 
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where we define 


K63(t) 

^53(t) 

MW 

v(t) 


D +(t) Jo 

7 

W(t') + 

( l D_{t)D + {t')-D + {t)D_(t') 

2 D 3 + (t) J 
21 

0 W(t>) 

D_{t)D+{t')-D + {t)D_(t') 

2 D\{t) J 

9 I'*' . 

0 W(t') 

D-(t)D + tf) - D+(t)D-(f) 3 


Dl(t')N aa (t')dt', 




(72) 

(73) 

(74) 

(75) 


Here note that the parameters in front of W a7 ij(p) and W a7 i(p) in expression (71) are the same, which 
originates from the relation (69). In the limit of the Einstein de Sitter universe in the general relativity, 
the coefficients, Kg 3 (t), A <53 (i), n(t), and v{t) reduce to 1. 

We can redefine these coefficients using the differential equations. Inserting the general form of the 
solution (71) into (61), we obtain the following differential equations for the coefficients 

«« 3 (f) + (6/ + 2 )k S3 (t) + (6 f 2 H 2 - 2 L) KS 3 {t) = N aa (t), (76) 

A ss(t) + (6/ + 2)A 43 (i) + (6 f 2 H 2 - 2L)A 43 (i) = - 7 -N ai {t ), (77) 

m + (6/ + 2)A(t) + (6 f 2 H 2 - 2L)n(t) = -yiV 7Q (t), (78) 

v(t ) + (6/ + 2 )v{t) + (6 f 2 H 2 - 2 L)v(t) = 9 N^t). (79) 

The homogeneous solution of all these equations is 1/D?(t) and D_(t)/D+(t). Therefore, the differential 
equations (76) to (79) consistently yield the inhomogeneous solutions (72) to (75), respectively. 

We next show that KS 3 {t) = 1 identically. Using the expression (63), we easily find that K 53 = 1 is the 
solution of (76). This means that the inhomogeneous solution (72) reduces to K 53 = 1. We can prove 
KS 3 = 1 directly from (72), using partial integral. 

Furthermore we can show that Xs 3 (t) = A (t) identically. We can rewrite Eq. (77), as follows, 


Xs 3 {t) + (4/ + 2 )H\s 3 (t) + ( 2f 2 H" — L)Xs3(t) + 2fH(\s3 ~ A) 

+ ( 4 fH 2 - L)(Xs 3 - A) = 7 -(fH 2 - Nj), (80) 

where we used (69) and (65). We can easily check that Xs 3 (t) and A(t) satisfies the same differential 
equation (see Ecj. (A28)), which leads to X$ 3 (t) = X(t). 

In summary, we have the expression equivalent to (71), 


J) =Dl 




Thus the third order solution of density contrast is characterized by A(t), fi(t), and u(t). Note that A (t) 
is defined to describe the second order solution, then n(t) and v{t) are the new coefficients which appear 
at the third order. Table I summarizes the parameters and the mode-coupling functions necessary to 
describe the second order solution and the third order solution. 

Recently, the authors of [44] investigated the third order solution of the density perturbations, in a 
similar way, but within a model of the general relativity. In their paper, 6 parameters are introduced to 
describe the third order density perturbations. Our results suggest that less number of parameters are 
only independent. 
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parameters 

mode-coupling functions 

5 2 {t,p) 

A CO 

VVce(p) 

e 2 (t,p) 

Ae(0 

W 7 (p) 

S3(t,p) 

A(0, MO, MO 

W<m(p), W a7 J?(p), Wc 7 i(p) 

e 3 {t,p) 

A(0, Ae(0, Me(0, v e{t) 

w 77 (p), w £ (p) 


TABLE I: Functions for the mode-couplings and parameters necessary to describe the second order solution and 
the third order solution. 


Inserting the solution (81) into Eq. (59), we find the solution for the velocity divergence 


where we define 


03(t,p) = -D%(t)f(w aa (. P) - 7f\d(t)W al R - yA(t)W a7 i(p) 


- ^e(t)W t ^(p) + ^fl(t)Wf(p)j, 

(82) 

fX 9 (t) = MO + 3 

(83) 

MO - *•(*) + 3 )0 

(84) 


Here note that A g(t) is the parameter to describe the second order solution, and /xg(t) and vg{t) are the 
new parameters which appear at the third order. 


In summary, we first introduced nine mode-coupling functions in the third order equations, (54) and 
(55) with (49). We find the three identities (56), (57) and (58). Then, only six mode-coupling functions are 
independent in the nine ones. This conclusion that the number of the linearly independent mode-coupling 
functions is six can be proved by using the generalized Wronskian. The coefficients in front of W aa R and 
WaaL in equation (54) are the same, which leads to the final third order solution (81) and (82) expressed 
in terms of the five mode-coupling functions. 


IV. POWER SPECTRUM 

The third order solution of the density perturbations enable one to compute the 1-loop (second order) 
power spectrum. The second order matter power spectrum was computed by many authors [27-33, 35- 
37, 44], in general relativity and modified gravity models. We find the expression for the 1-loop order the 


power spectra of density contrast and velocity divergence by 

<<S(t,k i)6(t,k 2 )) = (2 7 r) 3 5 (3) (k 1 +k 2 )P M (t,fc), (85) 

(tffrkOflftka)) = (27r) 3 5 (3) (ki + k 2 )(-/)P M (i,fc), (86) 

= (2TT) 3 S^(k 1 + k 2 )f 2 Pgg(t,k), (87) 

where we use the notation k = |ki|. Some details of their derivations are described in appendix C, we 
here show the results, 

Pss{t, k) = D 2 + (t)P L (k) + D%{t) (p 5 ( f } (i, k) + 2 p£ 3 \t, k)) , (88) 

Pse(t, k) = D 2 + (t)P h (k) + D 4 + (t) (P 5 ( f } (t, k) + 2Pj e 13) (t, k)) , (89) 

Pee(t, k) = D\{t)P^k) + D\(t) (p£ 2) (t, k ) + 2P e (33) (t, kj) , (90) 
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where D^_(t)Pi J (k) is the linear matter power spectrum, and we define 

p sf \*> k ) = g 8 ( 27r )2 J drP ^i rk ) J i dxP L {k(l + r 2 - 2rx) 1 / 2 ) 

((7 - 4A)r + 7z + 2(2A - 7)ra 2 ) 2 
(1 + r 2 — 2ra;) 2 

p if+ fc ) = 98(2n) 2 J drPh ^ J 1 dxP ^i k ( l + r2 - 2ra;) 1/2 ) 

((7 - 4A)r + 7x + 2(2A - 7)to 2 )((7 - 8A 9 )r + 7x + 2(4A@ - 7)ra 2 ) 

(1 + r 2 — 2rx) 2 

p eT\t, k ) = 98^)2 J drP ^ rk ) J i ^ p l(A:( 1 + r 2 - 2rx) 1/2 ) 

((7 - 8A@)r + 7x + 2(4A@ - 7)rx 2 ) 2 
(1 + r 2 — 2rx) 2 


(91) 


(92) 


(93) 


and 


2 pp3 \t,k) = 


k 3 


252(2tt) 2 


p l(A:) J drP L (rk) 


12fx— - 2(21 + 36A + 22/i) + 4(84 - 48A - ll/r)r 2 


- 6(21 - 12A - 2/i)r 4 + 4 (r 2 - l) 3 ((21 - 12A - 2/i)r 2 + 2/x) In 


r + 1 

It- — II 


(94) 


2+4 M) = 


k 3 


252(2tt) 2 


P L (fc) J drP L (rk) 


1 


6(/i + 3 /i@)— — 2(21 + 36A + 11/i + 33 fio) 


+2(168 - 96A - 11/i - 33/x*)r 2 - 6(21 - 12A - n - 3fi e )r 4 


H—j (r 2 - l) 3 ((21 - 12A - n - 3fig)r 2 + (/u + 3 /j, s )) In 


r + 1 

It- — II 


2+4 M) = 


k 3 


84(2tt) : 


; p h{k) J drP L (rk) 


12/r e — - 2(7 + 12A + 22/zg) + 4(28 - 16A - 11 n e )r 2 


-6(7 - 4A - 2 fj,g)r 4 + — (r 2 - l) 3 ((7 - 4A - 2/r e )r 2 + 2/x e ) In 


r + 1 

Ir-II 


(95) 


(96) 


The third order solutions of the density contrast and the velocity divergence are described by 6 parameters 
in Table I. The 1-loop power spectra are described by 4 parameters, and they do not depend on v(t) and 
ugit) (see Table II). In deriving the 1-loop power spectrum, we find that the relation, 


£(k,qi,—qi) = 0, (97) 

holds, which prevents the 1-loop power spectrum from depending on v(t) and vg(t). Details are described 
in appendix C and D. 



parameters 

Pss 

A(t), n{t) 

Pse 

A(t), A»(t), Afl(t), He(t) 

Pee 

A (t), A g(t), fj,g(t) 


TABLE II: Summary of the parameters to characterize the 1-loop order power spectra Pss, Pse, and Pee, respec¬ 
tively. 
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V. TRISPECTRUM 


Here we present the expression for the matter trispectrum in the real space, which is defined by 

{6(t, k 1 )6(t, k 2 )S(t, k 3 )S(t, k 4 )) c = (2 7 r) 3 ^ 3 )(k 1 + k 2 + k 3 + k 4 )T(t, k 4 , k 2 , k 3 , k 4 ). 

Using the solution up to the third order of the density perturbations, we find 

T(t , k 4 , k 2 , k 3 , k 4 ) = D+(t) ^Tii 22 (ki, k 2 , k 3 , k 4 ) + Tin 3 (ki, k 2 , k 3 , k 4 ) 


(98) 

(99) 


where we define 


Tn 22 (t, ki, k 2 , k 3 , k 4 ) = 4Pn(/ci)Pn(fc 2 ) [Pii(|k 4 + k 3 |)P 2 (f, k 4 , —k 4 — k 3 )P 2 (t, k 2 , k 4 + k 3 ) 

+Pn(|ki + k 4 |)P 2 (t,k 4 , —k 4 - k 4 )P 2 (f,k 2 ,ki + k 4 )] + 5 cyclic terms,(100) 
T m3 (f,ki,k 2 ,k 3 ,k 4 ) = 6Pii(/ci)Pii(fc 2 )Pn(fc 3 )P 3 (t, k 4 ,k 2 ,k 3 ) + 4 cyclic terms. (101) 

Note that Tn 22 contains P 2 , defined by Eq. (C7), which depends on A (f), while Tm 3 contains P 3 , defined 
by Eq. (C9), which depends on A (t), fi(t) and v{t). Therefore, the matter trispectrum depends on these 
three parameters. 


VI. APPLICATION OF KGB MODEL 


In this section, we exemplify the effect of the modified gravity on the 1-loop power spectrum. We here 
consider the kinetic gravity braiding (KGB) model [60, 66], which is considered in Ref. [26] to demonstrate 
the effect of the modified gravity on the bispectrum. We briefly review the model. The action of the KGB 
model is written as 


S = 



R + X) - G 3 (0, X)U^ + £ m 


( 102 ) 


where M p \ is the Planck mass, which is related with the gravitational constant Gjv by 87 tGn = l/M^. 
Comparing this action (102) with that of the most general second-order scalar-tensor theory (1), the 
action of the kinetic gravity braiding model is produced by setting 


G 4 = 


K 

2 


g 5 = 0, 


(103) 


and we choose K and G 3 as 


K = -A, 



(104) 


where n and r c are the model parameters. Useful expressions of the kinetic gravity braiding model are 
summarized in appendix E. 

When we consider the attractor solution, which satisfies i(j)HG^x = 1, the Friedmann equation is 
written in the form 


§-) =% + (l-fio) 
H 0 J a 6 


H 

Pn 


— 2/ (2n —1) 


(105) 


where Hq is the Hubble constant and Ho is the density parameter at the present time, and the model 
parameters must satisfy 


H 0 r c 


/ 2 n— 1 \ 

1 

V 3n ) 

[e(i — h 0 )J 


(2n—l)/4n 


(106) 
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On the attractor solution, L(t), defined by Eq. (32), reduces to 


L{t) 


3 2n + (3 n — l)S2 m ^. 2 
2 5 n — 


(107) 


where is defined by fl m (a) = flo-ffg/ H (a) 2 a 3 . The linear growth factor D + is obtained from Eq. (31) 
with (107) and (105). However, note that the quasi-static approximation on the scales of the large scale 
structures holds for n £ 10 (see [66]). 

The second order solution and the third order solution are obtained with 


Nj(t) = ii7 2 cr$(t) 


9 (1-On) (2 n-Sl m ) a 2 
4 D m (5n — fl m ) 3 


(108) 


H 2 ^(t) 


9(1 - On)(2 n - D m ) 3 (4n 2 (21 + 25AD m ) - 4nU m (21 + 10AD m ) + 12^(21 + 4AU m )) r2 

2812^(5 n-I2 m ) 5 


v$(t) = 0 . 


(109) 

( 110 ) 


We have A (t) from (A23) with (108). Using these results and Eqs. (66) and (67), we have the expressions 
for fi{t) and v(t) from (74) and (75). Eqs. (A24), (83), and (84) give expressions for A g{t), fig(t), vg(t), 
respectively. 

Table III lists the numerical values of these variables at the redshift z = 1, 0.5 and 0, for the KGB 
model with n = 1, 2, 5, as well as the ACDM model. 

Figure 1 shows A, fi, v, A g, fig, vg as function of the scale factor a. In each panel, the blue dash-dotted 
curve is the ACDM model, and the red dotted curve, the yellow dashed curve, and the green thick solid 
curve are the KGB model with n = 1, 2, and 5, respectively. All the curves take the limiting value 
unity at a = 0, but deviate from the unity as a evolves. Note that the deviation of A, fi, v from unity 
is small, of the order of a few percent, but the deviation of A g, vg is rather large, which could be 10 
percent. This is because the parameters associated with the velocity, A g and vg defined by Eqs. (A24) 
and (84), respectively, contain the time derivative term, which makes a large contribution. Plus, some 
part of the difference between the ACDM and the KGB model come from the difference of the growth 
rate /. Deviation of fig in the KGB model from that in the ACDM model is rather small compared with 
the deviations of A g and vg , which comes from the fact that fi is not a monotonic increasing function but 
there exists a maximum value at a ^ 1. 

Figure 2 shows the 1-loop power spectra Pss, Pse , Pee , from the top to the bottom, respectively, which 
are normalized by those of the ACDM model. These are the snapshots at z = 0, and we adopted the same 
normalization begin = 0.8 for each model, which means that all the models have the same linear matter 
power spectrum. In each panel, the red dotted curve, the yellow dashed curve, and the green thick curve 



ACDM 

KGB(n = 1) 

KGB(n = 2) 

KGB(n = 5) 

D+(z = 1 / 0.5 / 0) 

0.477 / 0.602 / 0.779 

0.496 / 0.642 / 0.858 

0.489 / 0.628 / 0.838 

0.484 / 0.620 / 0.827 

f(z = 1 / 0.5 / 0) 

0.869 / 0.749 / 0.513 

0.951 / 0.835 / 0.593 

0.919 / 0.813 / 0.605 

0.904 / 0.805 / 0.612 

X(z = 1 / 0.5 / 0) 

0.999 / 0.997 / 0.994 

1.000 / 0.999 / 1.003 

1.000 / 1.000 / 1.011 

1.000 / 1.002 / 1.019 

fj(z = 1 / 0.5 / 0) 

0.999 / 0.998 / 0.996 

1.000 / 1.001 / 1.015 

1.001 / 1.005 / 1.015 

1.003 / 1.007 / 1.011 

v(z = 1 / 0.5 / 0) 

0.998 / 0.996 / 0.991 

1.000 / 0.999 / 1.014 

1.000 / 1.003 / 1.034 

1.002 / 1.008 / 1.049 

Xg(z = 1 / 0.5 / 0) 

0.994 / 0.991 / 0.983 

0.998 / 0.995 / 1.043 

0.999 / 1.004 / 1.073 

1.003 / 1.014 / 1.095 

fie(z = 1 / 0.5 / 0) 

0.997 / 0.995 / 0.991 

1.000 / 1.006 / 1.041 

1.006 / 1.018 / 1.008 

1.010 / 1.021 / 0.974 

vg(z = 1 / 0.5 / 0) 

0.994 / 0.990 / 0.980 

0.998 / 0.998 / 1.089 

1.002 / 1.014 / 1.136 

1.009 / 1.030 / 1.169 


TABLE III: Numerical values of the growth factor D+, the linear growth rate f, and the coefficients A, fi, v, 
A q, fig, vg at the redshifts z = 1.0, 0.5 and 0, for the ACDM mode and the KGB model with n = 1, 2, 5. In each 
cell, a set of the three numerics means the values at the redshift 2 = 1.0, 0.5 and 0 from left to right, respectively. 
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A (a) 



^e(a) 



a 


M( a) 


Me( a) 




17 ( a ) \z 0 (a) 




FIG. 1: A, p, i/, Ag, /rg, i/g as function of the scale factor a. In each panel, the blue dash-dotted curve is the 
ACDM model, and the red dotted curve, the yellow dashed curve, and the green thick solid curve are the KGB 
model with n = 1, 2, and 5, respectively. 
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k(hMpc~') 



"00 KGB, 


r 00ACDM 



k(liMpc *) 


FIG. 2: Relative deviation of the power spectra Pss(k) (top panel), Pse(k) (middle pane.), Pee{k ) (bottom 
panel), under the kinetic gravity braiding model with n = 1 (red dotted curve), n = 2 (yellow dashed curve), 
n = 5 (green thick curve), which are divided by those under the ACDM model. The panels show the snapshot at 
the redshift z = 0. 
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show the KGB model with n = 1, 2 and 5, respectively. In the linear regime k ;$ 0.1[/iMpc _i ], all the 
models converge because they have the same linear matter power spectrum due to the same normalization 
as = 0.8. The differences between the KGB model and the ACDM model appear for the quasi-nonlinear 
regime k 0.1[/iMpc -1 ] due to the nonlinear effect. Because all the model have the same linear matter 
power spectrum, this figure shows that the enhancement of the power spectrum due to the nonlinear effect 
is small in the KGB model compared with that in the ACDM model. This is understood as the results 
of the Vainshtein effect. Furthermore, the deviation from the ACDM model is more significant in the 
velocity power spectrum than that in the density power spectrum. In general, the amplitude of the 1-loop 
power spectra Pgg, Pgg , and Pgg are decreased when any of A(f), //(f), A g(t), and /tg(f) is increased. The 
behavior of Pse and Pgg in the quasi-nonlinear regime is dominantly influenced by A g{t) and /J.g(t). 


VII. SUMMARY AND CONCLUSIONS 

We found the third order solutions of the cosmological density perturbations in the Horndeski’s most 
general scalar-tensor theory assuming that the Vainshtein mechanism is at work. We solved the equa¬ 
tions under the quasi-static approximation, and the solutions describe the quasi-nonlinear aspects of 
the cosmological density contrast and the velocity divergence under the Vainshtein mechanism. In this 
work, we thoroughly investigate the independence property of the mode-couplings functions describing 
the non-linear interactions. We found that the third order solution of the density contrast is characterized 
by 3 parameters for the nonlinear interactions, one of which is the same as that for the second order 
solutions. The third order solution of the velocity divergence is characterized by 4 parameters for the 
nonlinear interactions, two of which are the same parameters as those of the second order solutions. The 
nonlinear features of the perturbative solutions up to the third order are characterized by 6 parameters. 
Furthermore, the 1-loop order power spectra obtained with the third order solutions are described by 4 
parameters. Assuming the KGB model, we demonstrated the effect of the modified gravity in the 1-loop 
order power spectra at the quantitative level. We found that the deviation from the ACDM model ap¬ 
pears in the power spectra of the density contrast and the velocity divergence, which can be understood 
as the results of the Vainshtein mechanism. The deviation from the ACDM model is more significant in 
the velocity divergence than the density contrast, which is explained by a dominant contribution of the 
parameters A g(t) and fig(t). It is interesting to investigate whether this is a general feature of the modified 
gravity with the Vainshtein mechanism or not. 
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Appendix A: Coefficients in equations and solutions 


We summarize the definitions of the coefficients in the gravitational and scalar field equations (5), (6), 
(9). 


Ao 

Ax 

A 2 


Bo 


B 1 

B 2 

Bo 

Co 


e e g T s + v 

~ Jp + h +Tt ~ 2Gt ~ 2 1i-^2ip' 

= j [ + Gt-Bt, 

= g T - 

1 H' 


X 

H 


|</>G 3 a + 3 ^A + 2HX^j Gixx + 2AAG 4 axa - StfiG^x + 2<pXG 4 ^xx 


+ (H + H 2 '} cfrG 5x + ■ 


2 HX 


(h + A 2 ) 


X 


G-,xx + Hcj)XXGoxxx 


—2 (x + 2 Hx) G 50 a - ^IG 5 #x -x(x- 2 Hx) G 50 aa j, 


= 2A 


Gax + & {G 5X + XGoxx ) — G 50 + XGo<j,x 


= —2X 


[Gax + 2 A"G 4 xx + H(j)Gox + HcjiXGoxx — G 5^, — XG^x'j , 


= Hcj)XGoxi 


= 2X‘ 


oy2 / v 

G 4 AA H-g— \ 2<f>G5xx + (pXGoxxx — 2Gs^,a + XG^xxj 


Gi = Hcj)X (G 5 a + XGoxx ), 


(Al) 

(A2) 

(A3) 


(A4) 

(A5) 

(A6) 

(A7) 

(AS) 

(A9) 


where we define 


Jr 




2 

2 


G4 — A (jf>Gox + Gs^ , 

G 4 - 2AG 4 a - A (H<j>G 5X - G 50 )] , 


0 = —^AG 3 a + 2HG 4 - 8HXG iX - 8HX 2 G 4X x + <j>G H + 2X<j>G 4< /, x 
—H 2 <p (5AG 5 a + 2A 2 G 5 aa) + 2AA (3G 50 + 2XG^x ), 

5 = 2AA a - A + QX^HGox - 2AG 30 - 6A 2 G 4 + 24A 2 A(G 4 a + AG 4 aa) 
-\2HX4>G^x - 6A0G 40 + 2H 3 X<fi(5G 5X + 2AG 5 aa) 

-6A 2 A(3G 50 + 2AG 50 a), 


V = K — 2A(G 3 0 + <j)Gox) + 2(3A 2 + 2A)G 4 - 12A 2 AG 4 a - 4HAG 4 a 

-8HXG 4 x - 8HXXG 4X x + 2(4> + 2Hj>)G A4> + 4AG 4 ^ + 4A(<£ - 2H^)G A4>X 
—2X(2H 3 <j) + 2HH<j> + 3H 2 4>)G 5X - AH 2 X 2 4>G 5X x + 4iJA(A - HX)G 5<t>x 
+2 [2(HX)-+3H 2 X] G 50 + AHXjxGw 


(A10) 

(All) 

(A12) 


(A13) 


(A14) 












The coefficients in the first and the second order solutions are defined as follows, 
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m 

S(t) 

m 

Z(t) 

«q(*) 

Ny(t) 

m 


AqFt — Af, 

AqQt + A 1 A 2 , 

A\Qt + A 2 T t , 

2 (-Ao^t + 2AiA 2 Gt + A 2 J~t) , 

fhnTZ 

H 2 Z ’ 

/5m<^ 

H 2 Z 1 
PmT 
H 2 Z 1 
i/4 

-( 2B 0 Kq — 3 Bik^Kq — 3 B 2 k$Kq — 6 


D_{t)D + (t')-D + (t)D_(t') 


2 D\(t) 


W(t') 


D 2 + {t') (. f 2 H 2 -N 7 (t'))dt 


Ae(t) — A(f) + 
A$ (t) 

A«p (t) 


M 

2 fH' 


2 1 

—-K$A(i) + — ( 2BqTkq — 3BiSkq — 3B 2 1Zkq — QB^TZk^hq) , 


7 Z 

— KvpA(t) H—— (2BqA 2 GtHq + B\(A 2 kq 2A2t/TK>i/KQ) B 2 (Skq 


— 2B 3 (Sk< S ,Kq — A 2 K<pKQ + 


A<g(t) — ——KQA(t) — — (^2 BqG^kq ~t~ Bi(A 2 Gt k q ~ 2G^K\j,kq) ~h B 2 (7 ~kq — 
-j-2B 3 (7~Kq,KQ + A 2 Gt k & k q ~ Gt k ® k '&)) ■ 
are described in the previous paper [26], but one can show that A(t) 


Some details 
equation, 


A(t) + (4/ + 2)H\(t) + (2 f z H 2 - L)X(t) = -( f 2 H 2 - iV 7 ). 


(A15) 

(A16) 

(A17) 

(A18) 

(A19) 

(A20) 

(A21) 

(A22) 

(A23) 

(A24) 

(A25) 

+ 2A 2 Gtk®kq) 

(A26) 

2QrpK^KQ) 

(A27) 

obeys the differential 
(A28) 
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The coefficients for the third order solutions are defined as 
2 

cr<j>(t) = — (2BqTkq — 3B±Skq — 3 B 2 TIkq — QBzHk^kq) , (A29) 

H<b(t) = — (2B$Tkq\q — B\ (2Skq\q + TkqX^) — B 2 (21Zkq\q + T kqX<&) 

—25 3 (TZk^Xq + 'RsKqXxij + <SavqA$)) , (A30) 

^ (-3c 0 r4 - 4 Ci^4) . ( A31 ) 

2 

C\p(t) = — ( 2 B 0 A 2 GTkq ~ B\ (2A2Gtk<siKq — AJkq) — -B 2 (Skq + 2A2Gt k ®kq') 

— 2 U 3 (tSK\jy kq T A 2 GTk\£i A 2 K<p kq^J , (A32) 

2 

= — (2B 0 A 2 GtkqAq — B 1 (A 2 Gt (k^Aq + kq A$) — A^kqAq) 

—i ?2 ( SkqAq + A 2 GT (h<s>Aq + KQ A$)) 

B 3 (5 (k>i>Aq + kqAvp) + ^ 2 ^ 7 ’ (k$Avp + k>[/A$) A 2 (k$Aq + kqA$))) , (A33) 

^(t) = ( — 3CoA 2 t/TKQ — Ci (5 acq + 3A 2 C?tk<i>kq)) j (A34) 

CTQ(t) = ~ (— 2,BoGt^Q 4“ -®1 ( 207 ’KvpAtQ — A 2 GTKQ ) 1?2 (TKq 2Q^k^kq^ 

2 B 3 (Skq Qrp k^K- q/ T A 2 GTkqkq')') , (A35) 

2 

MQW = ^ (—2B 0 Gt k qAq + B\ (G? (k^Aq + kqA$) — A 2 GtkqAq) 

—B 2 {TkqAq — Gt ( K i>A q + K Q A$)) 

— B 3 (( SkqAq + TkqA$i) — Gt (k$A\J( + K\pAij>) + A 2 GT {k^Aq + kqA$))) , (A36) 
v q{^) — 7 ^ ( 3CqGt k q + Ci ( — ’ Tkq + 3 Gt k, s> k q)) ■ (A37) 

Appendix B: The third order mode-coupling functions 

In this appendix, we summarize the functions that describe the nonlinear mode-couplings of the third 
order solutions. In order to derive Eqs. (46), (47) and (48), we define 

W 7 a(p) = ^L_ / dkidk 2 dk 3 ^>(ki + k 2 + k 3 - P) 7 a(ki, k 2 , k 3 )^(ki)^(k 2 )^(k 3 ), (Bl) 

W 77 (p) = j tfcidkadM^k! +k2+k 3 -p)77(k 1 ,k 2 ,k 3 )^(k 1 )J L (k 2 )5 i (k 3 ), (B2) 

W € (p) = ydk 1 dk 2 dk 3 d( 3 )(k 1 +k2 + k 3 -p)e(k 1 ,k2,k 3 )d i (k 1 )d i (k 2 )<5 L (k 3 ), (B3) 


7 a(ki,k 2 ,k 3 ) = - ^ 7 (ki,k 2 + k 3 )a (s) (k 2 ,k 3 ) + 2 cyclic termsj , 

77 (ki, k 2 , k 3 ) = i ( 7 (k 1; k 2 + k 3 ) 7 (k 2 , k 3 ) + 2 cyclic terms), 

^ 1 ^i(k 2 • k 3 ) 2 + fc|(k 3 • ki) 2 + fc 3 (k! • k 2 ) 2 - 2(kx • k 2 )(k 2 • k 3 )(k 3 • kx) 

s(k;.k 2 .k :! i - I ---777-2- 

Aii " / 2 /i 3 


(B4) 

(B5) 

(B6) 


with 
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In deriving Eqs. (54) and (55), we define 

W aQfl (p) = J dkidk 2 dk 3 (5 (3) (ki+k 2 +k 3 -p)aai{(ki,k 2 ,k 3 )(5 L (ki)<5 L (k 2 )<5 L (k 3 ), (B7) 

W Q7fl (p) = y'dk 1 dk 2 dk 3 (5 (3) (k 1 +k 2 +k 3 -p)a7 fl (k 1 ,k 2 ,k 3 )(5 z ,(k 1 )(5 L (k 2 )(5 i (k 3 ), (B8) 

W qq l(p) = y'dk 1 dk 2 dk 3 (5 (3) (k 1 +k 2 +k 3 -p)aa i (k 1 ,k 2 ,k 3 )(I L (k 1 )5 L (k 2 )(5 L (k 3 ), (B9) 

W a 7 L( P ) = j dk 1 dk 2 dk 3< 5( 3 )(k 1 +k 2 +k 3 -p)a7 L (k 1 ,k 2 ,k 3 )5 L (k 1 )5 L (k 2 ) ( 5 i (k 3 ), (BIO) 

W QQ (p) = J rfk i dk 2 dk 3^ (3) ( k i+ k 2+k 3 -p)aa(ki,k 2 ,k 3 )5 L (ki)(5L(k 2 )(5 L (k 3 ), (Bll) 

W Q7 (p) = y'dk 1 dk 2 dk 3 (5 (3) (k 1 +k 2 +k 3 -p)a7(k 1 ,k 2 ,k 3 )(5 i (k 1 ) l 5 L (k 2 )5 L (k 3 ) (B12) 

with 


aa_R(ki,k 2 ,k 3 ) 

= ^ (a(ki,k 2 + k 3 )a (s) (k 2 ,k 3 ) + 2 cyclic terms) , 

(B13) 

a7_R(ki, k 2 , k 3 ) 

= i (a(ki, k 2 + k 3 ) 7 (k 2 ,k 3 ) + 2 cyclic terms), 

(B14) 

aa L ( ki,k 2 ,k 3 ) 

= i (a(ki + k 2 ,k 3 )a (s) (k 2 ,k 3 ) + 2 cyclic terms) , 

(B15) 

a7i(ki,k 2 ,k 3 ) 

= ^ (a(ki + k 2 , k 3 ) 7 (k 2 , k 3 ) + 2 cyclic terms), 

(B16) 

aa(ki,k 2 ,k 3 ) 

= i (a^(ki,k 2 + k 3 )a( s) (k 2 ,k 3 ) + 2 cyclic terms) , 

(B17) 

a7(ki,k 2 ,k 3 ) 

= i ^a (s) (ki,k 2 + k 3 ) 7 (k 2 ,k 3 ) + 2 cyclic terms) . 

(B18) 


Appendix C: Derivation of the 1-loop power spectra 

The cosmological density contrast S(t, k) and the velocity divergence 6(t, k) up to the third order of the 


perturbative expansion are expressed as 

S(t, k) = £>+(*)&(k) + D 2 + (t)6 2K (t, k) + D 3 + (t)6 3K (t, k), (Cl) 

0(t, k) = -/ (D+(t)S h ( k) + D 2 + (t)9 2K (t, k) + D 3 + {t)6 3K {t, k)) , (C2) 

were we define 

5 2K (t,k) = Wa(k) - ^A(i)W 7 (k), (C3) 

5 3K (t, k) = W« Q (k) - yA(<)W Q7 ij(k) - yA(t)W Q7 i(k) - ^-/r(t)W 77 (k) + ii/(i)Wg(k), (C4) 
6 2K (t, k) = Wa(k) - ^A 9 (i)W 7 (k), (C5) 


9 3K {t, k) = W aa (k) - yA e (t)W Q7fl (k) - yA(t)W Q7i (k) - y/r e (f)W 77 (k) + i^(t)W c (k), (C6) 
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and the kernels for the density contrast P 2 , and P 3 , and those for the velocity divergence G 2 , and G 3 , as 
follows, 

F 2 (t, ki,k 2 ) = a(ki,k 2 ) - ^A(i) 7 (ki,k 2 ), (C7) 

G 2 (t, ki,k 2 ) = a(ki,k 2 ) - ^Ae(t) 7 (ki, k 2 ), (C 8 ) 

2 2 

F 3 (t, ki,k 2 ,k 3 ) = aa(ki,k 2 ,k 3 ) - -A(i)a 7 jj(ki, k 2 ,k 3 ) - -A(i)a 7 i,(ki,k 2 ,k 3 ) 

- ^ A< ( t ) 77 ( k i’ k 2 ’k 3 ) + ^(i)C( k i, k 2 ,k 3 ), (C9) 

4 2 

G 3 (t,ki,k 2 ,k 3 ) = aa(ki,k 2 ,k 3 ) - -Ae(t)a 7 fl (ki, k 2 , k 3 ) - -A(t)a 7 z,(ki, k 2 , k 3 ) 

-^e(*) 77 (ki,k 2 ,k 3 ) + iv 0 (i)£(ki, k 2 , k 3 ). (CIO) 

These kernels have the two types of symmetries. One is the symmetries in replacement of the wave 
numbers, 


Flit, k]., k 2 ) — F 2 (t, k 2 , kx), (Cll) 

ki,k 2 ,k 3 ) = -F 3 (t, k 2 , k 3 , k 3 ) = F 3 (t, k 3 , k 1; k 2 ) 

= F 3 (t, k 1; k 3 , k 2 ) = P 3 (t, k 2 , k 3 , k 3 ) = P 3 (t, k 3 , k 2 , k 3 ). (C12) 


The second is the symmetries in the conversion of the sign of the wavenumbers, 

F 2 (t, ki,k 2 ) = F 2 (t, — ki, — k 2 ), (C13) 

Fi(t : k-|, k 2 , k 3 ) = F 3 (t, —ki, —k 2 , —k 3 ). (C14) 

The same relations hold for G 2 (t, ki,k 2 ) and G 3 (f, ki, k 2 , k 3 ). 

The above properties are useful in deriving the expressions of the power spectra, Pss(t,k), 
Pse{t,k), Pgg(t,k), defined by Eqs. (85), (86) and (87). Using the expressions (Cl) and (C2), we find 

Pss{t,k) = D 2 + (t)P h (k) + D%{t) (P^f\t,k) + 2P& 3 \t,kj) , (C15) 

Pse(t,k) = Dl(t)P L (k) + D 4 + (t)(p£ 2 \kk) + 2P£ 3 \t,k)) , (C16) 

Pee(t,k) = D 2 + (t)P L (k) + D 4 + (t)^ 22 \t,k) + 2P^ 3 \t,k)) , (C17) 

where is the linear matter power spectrum, 

<<5 L (k 1 ) < J L (k 2 )) = (2 7 r) 3 <5( 3 )(k 1 + k 2 )P L (fc), (C18) 

and we define 

(6 2 K(tM)S2K(t,k 2 )) = (27T) 3 ^(k 1 +k 2 )pjf ) (t,fc), (C19) 

(Mkr)^,^)) = (2 7 r)V 3 )(k 1 + k 2 )p£ 3) (i, fc), (C20) 

kr)^^(t, k 2 )) = (2 7 r) 3 ^ 3 )(k 1 +k 2 )p£ 2) (t,k), (C21) 

(e 2K (tM)0 2 idt,k 2 )) = (2n) 3 6^ 3 \k 1 +k 2 )P^ 2 \t,k) : (C22) 

(S^k^KitM)) = (2n) 3 5^(k 1 +k 2 )P^\t,k), (C23) 


2 ((<^L(ki )0 3 K{t, k 2 )) + (S 3 K(t, ki)^L(k 2 ))) 


(27T) 3 ^ 3 )(k 1+ k 2 )P] e 13) (t,fc). 


(C24) 


and 
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( 22 ) 

As an example, let us explain the derivation of Pg S '(t,k ). Inserting (C3) into (C19), we have 
{S 2 K{t,ki)S 2 K(t, k 2 )} = ^(2^)3 /AlA 2 ^( k l-qi-q 2 )^(,qi,q2).L(qi).L(q2) 

x (2^)3 J 94 ^ (3) ( k 2 -q3 - q4)^2(i,q3,q4)^L(qs)(5L(q4)^ 

= J d 3 q 1 d 3 q 3 F 2 (t, quki - q 1 )F 2 (t,q 3 ,k 2 - q 3 ) 

x (<5 L (qi)<5L(ki - qi)<5L(q 3 )^(k 2 - q 3 )} ■ (C25) 
Using the relation that hold for the Gaussian variables, we have 

(<5 L (qi)M k i - qi)<5i(q 3 )^L(k 2 - q 3 )} = (5z,(qi)^(ki - q 3 )} (<5 L (q 3 )<5L( k 2 - q 3 )) 

+ (^z,(qi)^(q3)} (^(k 2 - q 3 )<5 L (ki - qi)} 

+ <^(qi)^i(k 2 - q 3 )} (<5i(ki - qi)5i(q 3 )} , (C26) 

which yields 

(^z,(kn)<5L(ki — kn)5L(k 2 i)5i(k 2 — k 23 )} 

= (27r) 6 d (3) (qi + q 3 )6 (3) (ki + k 2 - qi - q 3 )P L (< 7 i)PL(|ki - qi|) 
-t-(27T) 6 A( 3 >( qi +k 2 - q 3 )<5 (3) (k! - qi +q 3 )P L ( gi )P L (g 3 ), (C27) 

with Eq. (C18). Then, (C25) yields 

(S 2K (t,k 1 )S 2K (t,k 2 )) =<5 (3) (k 1 +k 2 ) J d 3 q 1 (F 2 (t,q 1 ,k 1 - qi)P 2 (t, -qi, k 2 + qi) 

+P 2 (t,qi,ki - qi)P 2 (t,k 2 + qi,-qi)^P L (gi)P L (|ki — qi|). (C28) 

Using the relation (C13), we have 


p SS 2) ( t ^ k ) = J d 3 <?iP 2 2 (t,qi,k-qi)P L (gi)P L (|k-qi|), 


(C29) 


which reduces to (91). In the derivation, we define x = cos0, where 9 is the angle between ki and q 3 . 
Similarly, the expressions (92) and (93) are obtained for P$g 2 \t,k) and Pgg 2 \t,k). In the limit of the 
Einstein de Sitter universe withtin the general relativity, A (t) = X$(t) = /r(t) = /i g(t ) = 1, which gives the 
well-known expressions 


4 22) (*o = 

k 3 

^ drP\,{rk) 

j c!xPl(( 1 + r 2 — 2 rx) 1 ! 2 ) 

98(2t r) 2 j 

p£ 2 \k) = 

k 3 

f drP\,{rk) 

j dxPi,((l + r 2 — 2 rx) 1 / 2 ) 

98(2tt) 2 j 


r . 2'\2 


Hf\ k ) = 


k 3 


98(2t r) ; 


■ J drP\,{rk) J dxPh((l + r 2 — 2 rx) 1 ^ 2 ) 


(3 r + 7x — 10rx 2 )(—r + 7x — 6rx 2 ) 
(1 + r 2 — 2 rx) 2 

(—r + 7x — 6rx 2 ) 2 


(1 + r 2 — 2 rx) 5 


(C30) 

(C31) 

(C32) 


which are constant as functions of time. 

Next, let us explain the derivation of Pg S (t,k ). Inserting (C4) into (C20), we have 

2 {SL(ki)6 3 K (t, k 2 )) 

2 


(2t r) 6 


J d 3 q 1 d 3 q 2 F 3 (t, qi,q 2 ,k 2 - q 3 - q 2 ) (<5L(ki)5 L (qi)<5 L (q 2 )5L(k 2 - qi - q 2 )). (C33) 
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Using the relations, 

(Mki)£i(qi)<Sz,(q2)£L(k2 - qi - q 2 )} 


and the symmetries, (Cl 2), we have 


2 P& 3 \t,k) 


6 


( 2 ^) 6 ^(k 1 + q 1 )P L (fc 1 )^ (3) (k 2 - qi)P L (g 2 ) 

(27r) 6 ^ 3) (ki +q 2 )Ph(fc 1 )^ 3) (k 2 -q 2 )i\( gi ) 

(27r) 6 <5 (3) (ki + k 2 - qi - q 2 )P L (ki)5 (3) (q : + q 2 )P L (o , i), 

(C34) 


d 3 qiF 3 (t, k, qi, -qi)P L (fc)P L (gi). (C35) 


After performing the angular integration with respect to the spherical coordinate of qi, we finally have 
(94). Note that (94) does not depend on v(t), which occurs because of the identity £(k, qi,— qi) = 0. 
Pgg 3 \t, k) is characterized by A (t) and Similarly, we have the expressions (95) and (96) for P^ 3 \t, k) 
and Pgl 3 \t, k), respectively. Because of the same reason for P^ 3 \t, k), Pgl 3 \t, k) and P g l g 3 \t , k) do not 
depend on v{t) and ve(t). Furthermore, P^ 3 \t, k) and P g l g 3 \t, k) do not depend on A g(t). This is because 
of the nature of the integration 


J dx cr^k, qi, —qi) = 0. 


(C36) 


Finally, P^ 3 \t,k) depends on A(f), /r(f), and ne(t), and P g l g 3 \t,k ) depends on A(t) and /rg(f). We find 
the following relation holds, in general, P gg 3 \t, k) = [P gg 3 \t, k) + P g l g 3 \t , k)\/ 2, from (C24). 

In the limit of the Einstein de Sitter universe withtin the general relativity, all the coefficients A(f), 
/r(t), /xg(t) reduce to 1, which reproduces the well-known expressions 


2 PjfW = 


2Pj e 13) (fc) = 


2P^ 3 \k) = 


k 3 


252(2tt) 2 
x 


Ph(k) J drP L (rk) 


- 158 + 100r 2 - 42r 4 + dL(r 2 - l) 3 (7r 2 + 2) In ' T + 1 


r - 1 


k 3 


252(2 

x 24^ — 202 + 56r 2 — 30r 4 + 2- (r 2 — l) 3 (5r 2 + 4) In ( r | ^ 
r 2 r 6 K v ' \ \r — 1| / 


k 3 


84(2tt) 2 
x 


P L (fc) J drP L (rk) 


12-^ - 82 + 4r 2 - 6r 4 + (r 2 - l) 3 (r 2 + 2) In ( r + j 

r i r i \r-l 


(C37) 


(C38) 


(C39) 


Appendix D: The integrations of mode-coupling functions 


Here we summarize the useful expressions, which are useful in deriving the 1-loop order power spectra, 

(r + x — 2 rx 2 ) 2 


cU(qi,ki - qi) = 
a(qi,ki-qi) 7 (qi,ki - qi) = 
7 2 (qi,ki - qi) = 


4r 2 (l + r 2 — 2 rx) 2 ’ 

(r + x — 2rx 2 )(—1 + a; 2 ) 
2r(l + r 2 — 2 rx) 2 
(-1+x 2 ) 2 


(1 -F r- 2 — 2 rx) 2 ’ 


(Dl) 

(D2) 

(D3) 
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and 


/ 2irk 3 f 

d 3 qi P L (rk)aa( k,qi,-qi) = / drP h (rk) 

j d 3 qi Pl (rk)a"/R (k, qi, —qi) =0, 

/ 2irk 3 I" 

d 3 qi P L (rfc)a 7 L (k,q 1 ,-q 1 ) = J drP L (rk) 

/ 2 t T k 3 r 

d 3 qi P L (rfc) 77 (ki,qi,-qi) = -jA- J drP L (rk) 

J d 3 qi P L (r/c)£(ki,qi,-qi) = 0. 


—2 + 16r 2 — 6r 4 + -^-(r 2 — l) 3 In (-——- 
r 6 \ r-- 1 


6 + 16r 2 - 6r 4 + hr 2 - l) 3 In ( 

r J \ r-1 


(D4) 

(D5) 

(D6) 


-6-^ + 22 + 22r 2 - 6 r 4 + A(r 2 - l) 4 In ( 

r 6 \ r-1 


(D7) 

(D8) 


Appendix E: Coefficients and in the KGB model 

In the KGB model, we find the coefficients in basic equations, 

Ft = M 2 , g T = M 2 , 

e = -nM p i(-J r ] j>X n + HM 2 h 


Ph 


0 = —n(2n + 1 )M pl ^ + 7/M 2 , 

£ = -A + 6nM p i ( ^ J <])HX n - 3i/ 2 M 2 , 


^ P 2 i 


V = -X~ 2nM p i (t)X n + (3M 2 + 2 H)M 2 h 


^2 = Bn = n—M, 


P 1 m2 


H \ M, 


A" 


pi 


Ai — Bi — B 2 — B3 — C 0 — Ci — 0, 


(El) 

(E2) 

(E3) 

(E4) 

(E5) 

(E6) 

(E7) 

(E8) 


L(t) = - 
A 7 (t) = 

= 


AoFrPni 

2(AqGt + A\Ft) ' 
B 0 A 3 P 3 p 2 m 
4 {A 0 g 2 + A 2 F t ) 3 H 2 ’ 
8B 0 T 3 p 2 m 8 B 2 g 2 T A p 3 m 

7 H 2 Z 3 H 4 Z 5 


(E9) 

(E10) 

(Ell) 


and the non-trivial expressions, 
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We use the attractor solution which satisfies 2>4>HGox = 1- Then we have 


* = - 1 

Y 2n — 1 H 

(E12) 

H (2 n - l)3fi m 

H 2 “ 2(2 n - O m ) ’ 

(E13) 

M p 2 j( 1 - On) (2 n + (3 n - l)fi m ) 

0 “ 2n-n m 

(E14) 

A 2 = Mpj(l - Sl m ), 

(E15) 

B 0 = M 2 ( 1 - Q m ), 

(E16) 

where we define = p m {a) /3M pl iL 2 . We also have 


M£( 1 - On) (2 n + (3 n - l)fi m ) 

2 n — 

(E17) 

1 - n m ) (2 n + (3 n - l)fi m ) 

2 n — 

(E18) 

r = m p 4 !(i - n m ), 

(E19) 

^ o M® i n m (5n - n m )(l - fi m ) 

2 n — 

(E20) 
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